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Abstract 


A combined finite element method ( FEM ) and method of moments 
( MoM ) technique is presented to analyze the radiation characteristics 
of a cavity-fed aperture in three dimensions. Generalized feed model- 
ing has been done using the modal expansion of fields in the feed struc- 
ture. Numerical results for some feeding structures such as a 
rectangular waveguide, circular waveguide, and coaxial line are pre- 
sented. The method also uses the geometrical theory of diffraction 
(GTD) to predict the effect of a finite ground plane on radiation char- 
acteristics. Input admittance calculations for open radiating structures 
such as a rectangular waveguide, a circular waveguide, and a coaxial 
line are presented. Numerical data for a coaxial-fed cavity with finite 
ground plane are verified with experimental data. 

1. Introduction 

Cavity-backed aperture antennas are very popular in aerospace applications because of their confor- 
mal nature. The cavities can also have various types of microstrip antennas residing in them. Hybrid 
techniques have recently become attractive for the numerical analysis of these types of problems. The 
combined finite element method and method of moments (FEM and MoM) techniques, in particular, 
have been used to analyze various cavity-backed antenna systems (ref. 1). 

The FEM has a relatively simple formulation and is attractive for complex penetrable structures. 
Also, it results in sparse, banded matrices that can be effectively stored and solved. However, the FEM 
when used alone does not incorporate the Sommerfreld radiation condition and hence requires discreti- 
zation to extend far from the source region so that the radiation condition can be imposed. Recent 
efforts have concentrated on the use of absorbing boundary conditions to reduce the discretization 
region (ref. 2). Unfortunately, the accuracy of these approximate boundary conditions depends on spe- 
cific problem geometry, leading to results of unpredictable accuracy. 

The method of moments (MoM), on the other hand, incorporates the Sommerfield radiation condi- 
tion through the use of the appropriate Green’s function. As a result, the domain discretization can be 
kept to a minimum. However, this method has the disadvantage of being difficult to implement for com- 
plex penetrable structures. This method also results in full matrices whose treatment requires excessive 
storage, limiting its application to small-scale devices. 

To eliminate the disadvantages of both methods, several authors have proposed a combined 
approach (refs. 3 and 4) for radiation and scattering problems. For a comprehensive list of references on 
hybrid FEM/MoM techniques, see reference 5, chapter 9. 

The modeling of the coaxial feed in reference 1 is accomplished by replacing the feed probe by a 
constant filament of current with a delta-gap generator. The feed modeling in reference 1 is expected to 
give good results for shallow cavities, but deep cavities require rigorous feed modeling. In this paper a 
combined FEM/MoM approach is used to predict characteristics of cavity-backed antennas (fig. 1). 
Various feed lines such as a rectangular waveguide, a circular waveguide, and a coaxial line are consid- 
ered. Using the modal expansions of fields in the feed, the input admittance of cavity-backed antennas is 
calculated accurately. 

In most antenna-pattern measurements, the infinite ground plane is simulated by a relatively hu ge 
finite ground plane. Although this simulation does not appreciably affect the input admittance measure- 
ments, diffractions from the edges of the ground plane modify the radiation pattern considerably. In the 
past, the geometrical theory of diffraction (GTD) was successfully applied to predict the changes in the 
radiation pattern of the aperture antennas caused by a finite ground plane (ref. 6). The technique 



described in reference 6 predicts the changes in the radiation patterns of cavity-backed aperture anten- 
nas with finite ground planes. 

The body of this paper is organized in the following manner. Section 2 presents the theory for cav- 
ity-backed apertures in infinite ground planes. Section 3 compares the numerical results calculated 
using the present analysis with the available results in the literature. The input admittance calculations 
for an open rectangular waveguide, a circular waveguide, and a coaxial line are presented and compared 
with the results available in the literature. To verify the validity of the GTD technique used in combina- 
tion with the FEM/MoM technique, the radiation patterns of a coaxial cavity with finite ground plane 
are calculated and compared with experimental results measured at the Langley Low-Frequency 
Antenna Test Facility. Numerical data for input return loss of such a coaxial cavity are also computed 
and compared with experimental data. Section 4 presents concluding remarks. For completeness, an 
appendix includes expressions for modal expansions of electric fields of input-feed-line structures. 

Symbols 

A FEM/MoM matrix 

A p( j. B pq amplitude of modal vectors in feed 

a, b length and width of rectangular waveguide, respectively 

a m reflection coefficients for TEM/TE modes 

b excitation matrix 

b m reflection coefficients for TM modes 

D s , Dj ] diffraction coefficients 

ds surface integral 

dv volume integral 

E time-harmonic electric field 

E inc incident electric field 

E inp electric field coming from transmission line input 

e unknown coefficient matrix 

e t unknown coefficients as defined in equation (8) 

ej nc modal vector function for incident electric field in feed structure 

TEM/TE modes in feed line 

TM modes in feed line 
electric vector potential 
H time-harmonic magnetic field 

H ((/ , magnetic field at aperture 

H (/ total diffracted magnetic field 

diffracted magnetic field from point A 

Hg diffracted magnetic field from point B 

H^. diffracted magnetic field from point C 
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diffracted magnetic field from point D 
% far field of magnetic field 

h'u total magnetic far field after adding diffracted field 

Hj np magnetic field at input surface 

Hqu magnetic far field, 0-component 

H 0 jj magnetic far field, (^-component 

Hq a diffracted magnetic field from point A, 0-component 

A diffracted magnetic field from point A, (^-component 

Hq b diffracted magnetic field from point B, 0-component 

B diffracted magnetic field from point B, (^-component 

7| ( ) Bessel function of first kind (first order) 

J n ( ) Bessel function of first kind (nth order) 

J' ( ) derivative of J n ( ) with respect to its argument 

j = 'J-l 

k Q free-space wave number, ™ 

k 

kp cutoff wave number for TM modes in circular waveguide 

k'p cutoff wave number for TE modes in circular waveguide 

L distance between radiating aperture (in z = 0 plane) and feed-line input surface 

(in z-~L plane) 

Lj length of edge of tetrahedral element 

M magnetic current 

n wedge factor 

h unit normal 

/;, q mode indices for feed 

r position vector from origin to field point 

r' position vector from origin to source point 

r c radius of circular waveguide 

r 1? /'2 inner and outer radius of coaxial line, respectively 

r, 0, c() spherical coordinates 

S surface area 

S ap surface area of aperture (fig. 1) 

S inp surface area of input plane (fig. 1) 

T vector testing function for volume elements 

Tj vector testing function for surface elements 
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TE 

TE 10 

TE n 

TEM 

TM 
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a ct 
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P cA 
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P cC 
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p. (|), Z 
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transverse elective 

dominant mode in rectangular waveguide 
dominant mode in circular waveguide 
transverse electric and magnetic 
transverse magnetic 
volume 

vector basis functions for tetrahedral elements, i - 1, 2, . . . , 6 

vector basis functions for tetrahedral elements, j = 1, 2, ..., 6 

length and width of finite ground plane, respectively (fig. 1) 

rectangular coordinate axes (fig. 1) 

unit vector along X-, Y-, and Z-axis, respectively (fig. 1) 

input admittance of feed line 

characteristic admittance of feed transmission line 

simplex coordinates associated with nodes m and n in tetrahedral element 

angle between edge of diffraction and ray of incidence on edge 

reflection coefficient as defined in equation (27) 

propagation constant of incident field in feed 

propagation constant of feed line 

Neumann’s number 

relative permittivity 

intrinsic impedance of free space 

unit vector along 0-direction 

free-space wavelength 

free-space permeability 

relative permeability 

angle between ground plane and diffracted ray from point A 
angle between ground plane and diffracted ray from point B 
distance between center of aperture and diffraction point A 
distance between center of aperture and diffraction point B 
distance between center of aperture and diffraction point C 
distance between center of aperture and diffraction point D 
cylindrical coordinates 

unit vector along p-, (])-, and z-direction, respectively 
unit vector along (^-direction 
pth zero (p = 1, 2, 3 ...) of J n ( ) 
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l' pq Pth zero (p = 1 , 2, 3 . . .) of ( ) 

\|/ 0 angle between aperture (source) and point of diffraction 

to angular frequency 

A integration over a tetrahedral element 

V del operator 

2. Theory 

Figure 1 shows the geometry of the problem to be analyzed. For a linear, isotropic, and source-free 
region the electric field satisfies the vector wave equation: 

Vxf— VxEl-^E = 0 (1) 

vft f / 

where p,. and £,. are the relative permeability and relative permittivity of the medium in the cavity. The 
current formulation can be easily modified to incorporate anisotropic materials by using proper tensor 
quantities for u r and £,. (ref. 7). The time variation exp (/(Of) is assumed and suppressed throughout this 
paper. 

To facilitate the suitable solution of the partial differential equation in equation (1) with the FEM, 
equation (1) is multiplied by a vector testing function T, and the result is integrated over the volume of 
the cavity. By applying suitable vector identities, we can write equation (1) in its weak form (ref. 8) as 

UJi(V X T) • (V X EMv- de.JJjT E* = |||v • (t X lv X E )iv (2) 

V 1 V V ’ 

Applying the divergence theorem to the integral on the right-hand side of equation (2) and using 
Maxwell’s equation V x E = -joop 0 p 7 .H, where \i a is the permeability of free space and to is the angular 
frequency, we can write equation (2) as 

JJJ— ( VxT ) - ( VxE )z/v — ^eJJJt-EJv = jcopJJlTxn ) -Hds (3) 

V V s 

where h is the outward unit normal to the surface and H is the magnetic field at the surface. Because the 
tangential electric field is zero on the perfect-electric-conducting (PEC) walls of the cavity, the surface 
integral in equation (3) is nonzero only over the aperture opening in the infinite ground plane and the 
input aperture. Flence, we can write equation (3) as 

JJjl(VxT).(VxE)*- k 2 0 e r JJ|T • E dv - ,/copJJ (T x n) • H ap ds 

V ^ V S ap 

= jwpjj (Txn)-H iV A (4) 

^inp 

where S ap is the aperture surface and 5j np is the input surface (see fig. 1), H ((/ , is the magnetic field at the 
aperture, and H inp is the magnetic field at the input surface. 

At this point, the problem can be divided into three parts, the first paid involving the discretization 
and evaluation of volume integrals on the left-hand side of equation (4) and the second part involving 
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the discretization and evaluation of the surface integral over S ap . The third part involves the surface 
integral over the input aperture. These three parts are given explicitly as 


I v = jjj-J-(VxT)(VxE)dv-^8 r |||TEdv 
V V 

(5) 

1 ap = -^ojj (TX «)- V 

S ap 

(6) 

7 inp = Mkjj ( Tx «) ■ H 1|lf/ /x 

(7) 


^inp 


2.1. FEM Formulation 

To discretize the volume integrals in equation (5), the volume of the cavity is subdivided into small- 
volume tetrahedral elements. The electric field is expressed in terms of the edge vector basis functions 
(ref. 5), which enforce the divergenceless condition of the electric field implicitly: 

6 

E = Y. e i w / 181 

i = 1 

where e ,■ are the unknown coefficients associated with each edge of the tetrahedral element and 

W, = L ; .(a„ ; Va„-a ;1 Va m ) (9) 

where m and n are the nodes connected to form edge i, L ( - is the length of the z'th edge, and a m and a n are 
the simplex coordinates associated with nodes m and n (ref. 8). 

The testing function T is taken to be the same set of basis functions as that given in equation (8); 
that is, 


T = W j O’ =1,2,..., 6) (10) 

Substituting equations (8) and (10) into equation (5) and integrating over the volume of each ele- 
ment yields 

6 

t ev = -£-£ JJJ(VxW ; .) (VxW.) e .dv-^eJ||(W r W ; .) ei .dv (7 =1,2 6) (11) 

^ r i = 1 A A 

where ▲ indicates the integration over the volume of a tetrahedral element and I ev represents the inte- 
gral I v over a single element volume. 

The evaluation of these integrals over a tetrahedral element is given in reference 9. These element 
matrices are assembled over all the elements in the cavity volume to form a symmetric, sparse matrix of 
FEM. 

2.2. MoM Formulation 

The discretization of cavity volume into tetrahedral elements automatically results in the discretiza- 
tion of the aperture into triangular elements at the z = 0 plane. In accordance with the equivalence 
principle (ref. 10), the fields in the two regions can be decoupled by closing the aperture with a perfect 
electric conductor and introducing the equivalent magnetic current 

M = E x z (z = 0) (12) 
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over the extent of the aperture. Hence, from equation (8), the magnetic current over a triangular element 
can be expressed in terms of the unknown coefficients associated with the tetrahedral elements as 

3 

M = Y e,(W,xz)| (13) 

L -‘ lz = 0 

i = 1 

Similarly the testing functions can be modified as 

T = Txz = W-Xz (j =1,2, 3) (14) 

A J 

The magnetic field caused by the magnetic current given in equation (13) above the infinite ground 
plane is given by 


H fl = — (VxVxF) (15) 

ap i®p o 


where the electric vector potential F (making use of the image theory) is given by 


F 



exp(-jkjr-r'l) 


dx'dy' 


(16) 


where r is the field point and r' is the source point. 

Using the appropriate vector identities, we can rewrite the equation for the radiated magnetic field 
in equation (15) as (ref. 10) 


= ^7T [ ^ f + v ( v ' f )] d ? ) 

1 7®p o 0 

Substituting equation (17) into equation (6) and using standard surface vector calculus formulas 
(ref. 11) yields 

T ap = ~K\\ V F ^ + jj (V-T^V-F )ds (18) 

S a p $ap 

The integrals in the above equation are evaluated over all the triangles on the aperture surface fol- 
lowing the procedure outlined in Rao et al. (ref. 11). When calculating the mutual term over the trian- 
gles, a 13-point Gaussian quadrature formula developed for triangles (ref. 12) is used. Closed-form 
expressions given by Wilton et al. (ref. 13) are used to evaluate the self term of these integrals. The 
dense matrix thus formed is assembled over all the triangles and combined with the FEM matrix using a 
global numbering system. 

2.3. Input Integral 

For any transmission line input, the electric field can be formulated as the sum of the incident and 
reflected fields. 


w-, Tf 1 TM 

E inp = E inc + Y, [a m e m ( x >y) + b m e m U }’)]exp(Y m z) 
m 


(19) 


E inc = e inc eX P(-Tinc^) 


( 20 ) 
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TE TM 1 

where E inc is the incident electric field with unit amplitude, e m and e represent the TEM/TE and 
the TM mode, respectively, a m and b m are the reflection coefficients for the TEM/TE and the TM mode, 
respectively, and y m is the propagation constant. With the orthogonality properties of the waveguide 
modes (ref. 10), the reflection coefficients are 


a m = cx P(-Y m -)jj • ( E - E inc)^ d V 

^inp 

b m = exp(-y m z)JJ e™-(E-E inc )dx dy 

^inp 

The magnetic field at the waveguide input is given by 


H|np i®pip,. VXEin P 


( 21 ) 


( 22 ) 


(23) 


Substituting a m and b m in equation (23) and then using equation (23) in equation (7) and integrating 
over the feed-line cross section at z - ~L, we can write equation (7) as 


7 inp = 2 7,nc eX P(7,nc L )j{ T ' e inc dx 

^inp 

_ 1| T ' ^n dX J - V T || E ' e >n dX dy 


m V S 


inp 


inp 


II T ■ ™ dx dy II E ■ ™ dx dy 


m 1 m \ s 


(24) 


2.4. FEM/MoM Matrix Equation 

Combining all the integrals given in equations (5)-(7), we can write the system equation for this 
problem as 


JJ J — ( V x T) • ( V x E )dv - *o e /||| T ' E dv - /topJJ (T x h) ■ H ap ds 

T 7 T 7 C 


n 


T • zl E dx 


m V 5 


dy 11 : • e m dx dy 


inp 


inp 




m s 


inp 


inp 


= 2 Tinc eX P^inc L )|| T ' C mc dx dy 


(25) 


inp 


'Because the analysis is assumed to be valid for a generalized feed line, both TEM and TE modes are represented. They need 
only be used when applicable. 


8 



Substituting the expansion functions for E and T converts this equation into a set of linear equations 
to solve for the unknown cq given in equation (8). The linear equation system to be solved is written in 
matrix form as 


Ae = b (26) 

where A is a partly sparse (from the FEM) and partly dense (from the MoM and the modal expansion 
method) matrix, b is the excitation vector, and e is the unknown vector to be solved. The biconjugate 
gradient algorithm is used with diagonal preconditioning (ref. 14) to solve equation (26) efficiently. 
This algorithm also makes use of the symmetry of the matrix to store the matrix effectively. 

2.5. Reflection Coefficient and Input Admittance 

Solving equation (25) and using equations (21) and (22), we can obtain the reflection coefficients of 
any waveguide mode: 


r = a m or r = b m 


The input admittance is then calculated as 


Y = - — -y 

in 1 + p o 


where Y a is the characteristic admittance of the feed transmission line. 


(27) 


(28) 


2.6. Far Field With Infinite Ground Plane 

Once the electric field in the aperture is known, using equation (13) for the magnetic current, the 
magnetic far field with infinite ground plane can be evaluated directly from the following equation: 


H ff(r, 6 , 4 >) | 


r — > oo 


jk 0 exp(-jk 0 r) 

2n~r 


r r /v /v /v 

jj (00 + <M>) 


s 


ap 


■ M(x, y)exp [jk sin 0(x cos (|) + y sin (|))]dx dy 


(29) 


2.7. Application of GTD for Finite Ground Planes 

The infinite-ground-plane solution described in the preceding sections yields the field radiated in 
the forward direction. This solution is obtained with equation (29). The effect of the finite ground plane 
is incorporated by using the compact edge diffraction coefficients (ref. 15) and following the procedure 
described in reference 6. The radiation fields are supplemented by the additional field contributions 
arising from the diffraction by the edges of the finite ground plane. The diffracted field supplements the 
field in the forward region, whereas it accounts totally for the field in the shadow region (behind the 
finite ground plane) where the infinite-ground-plane solution gives a zero field. 

The geometry of the finite rectangular ground plane is illustrated in figure 1 . The midpoints on the 
two edges of length W3 are given by A and B, respectively. The midpoints on the two edges of length W\ 
are given by C and D . respectively. The diffracted fields from A and B are given by (ref. 6) 

< a - \n,Jd. 5 ojzyp^, q, v „ 

2 W e//(^’ 2 ' K A' - Vo- It, ■ " 1 


P cB 

^+P cA) 




exp (-jk r) exp(-jk d sin 0 ) 


(31) 


P cA 

r(r + Pc A ) 


-11/ ^ 


exp(-jk r)exp(jk dsin0) 


(30) 
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for cf> polarization, and 


< A - }*4 A ? °)°* (p ^ ti - K ■ "> 

2' % y >h ^ cB ’ ^ 2 ‘ 


PcA 

r(r + p cA ) 


-, 1/2 


exp(-j& o r)expO& o dsin0) 


P<~fi 

r ( r+ PcA) 


-] 1/2 


exp {-jk Q r) exp ( J sin 0 ) 


for 0 polarization. In the preceding equations, 


IT, 

PcA - P(« _ d - ~ 


(32) 


(33) 


(34) 


and 


= 5 + 0 for 0 < 0 < 7t 


^2 


5-0 for O<0<5 

2 2 


(35) 


(36) 


where 

\|/ 0 angle of incidence from aperture (source) toward the point of diffraction (A or R) = 0° 

71 

P' - for normal incidence (ref. 15, p. 811) 

n wedge factor (ref. 15, figs. 13-30) 

Similarly, the diffracted fields from C and D are also calculated by replacing p cC = p cD = d = W^. 
The detailed expressions for diffraction coefficients D s and D h are given by equation (113) of refer- 
ence 15. Total diffracted field is given by 


H 


d 


for «|) = 0° plane 
H d c + H d D for <> = 90° plane 


(37) 


The diffracted fields for the cf) = 180° and (f) = 270° planes can also be calculated in a similar man- 
ner. The total radiated far field is given by adding the diffracted field and magnetic far field given in 
equation (29); that is, 

H l jj = H ff + H d (38) 

The preceding analysis can be easily extended to a circular finite ground plane. If the circular 
ground plane is symmetrically excited, the diffracted rays come to a focus in the forward and the rear 
axial direction; that is, the axis of symmetry is a caustic of diffracted rays. The GTD field is not valid at 
the caustic. However, GTD can be used indirectly to calculate the equivalent edge currents (around the 
rim of the ground plane) that radiate in the caustic direction; these edge currents are integrated around 
the circular rim to yield the field in the neighborhood of the caustic. The radiated field evaluation is the 
same as that described in reference 6. 
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3. Numerical and Experimental Results 

To validate the present analysis, various numerical examples are considered. As a simple case, 
open-ended waveguide structures excited by a dominant mode field and radiating into half space 
through an infinite ground plane are analyzed. To study the effect of a finite ground plane, this paper 
considers a coaxially fed cavity that excites a coaxial aperture in a finite ground plane. The numerical 
data for these cases are compared with earlier published results and the experimental data obtained at 
NASA Langley Research Center. 

3.1. Apertures in Infinite Ground Plane 

3.1.1. Rectangular aperture. A rectangular aperture as shown in figure 2(a) with dimensions 
a = 4.755 cm, b = 2.215 cm, and L- 1.9 cm in an infinite ground plane and excited by a rectangular 
waveguide loaded with a dielectric plug (£,. = 3.76) is considered. The dielectric plug is discretized 
using 1230 tetrahedral elements. The air-dielectric interface inside the waveguide is considered as the 
input plane. The dominant TE 10 mode is considered to be incident at the plane z = —L. The input admit- 
tance of this structure is calculated as a function of frequency and plotted in figure 2(b). The numerical 
results agree very well with the data presented by Swift (ref. 16). 

3.1.2. Circular aperture. A circular aperture in an infinite ground plane and excited by a circular 
waveguide is considered as a second example. The geometry is shown in figure 3(a). The circular 
waveguide of length L = A from the aperture is considered as a circular cavity for FEM discretization. 
The waveguide is discretized at a level of 0.1 A,. The dominant TEjj mode is considered to be incident at 
the plane z = —L. The input admittance is computed as a function of the radius of the waveguide and 
plotted in figure 3(b) along with the numerical data given by Bailey and Swift (ref. 17). The results 
computed using the present technique and those of reference 17 agree very well. 

3.1.3. Coaxial aperture. As a third example, a coaxial aperture excited by a coaxial line is consid- 
ered. The geometry is shown in figure 4(a). Again the coaxial line of length L = A from the aperture is 
considered as a coaxial cavity for FEM discretization. The coaxial line is discretized at a level of 0.1 A. 
Only the dominant TEM mode is considered for both incident and reflected fields in the coaxial feed 
line because the higher order modes do not seem to contribute to the final result. (See appendix 
section A.3.) The input admittance is computed as a function of frequency and plotted in figure 4(b). 
The results are compared with the available data in the literature (ref. 18). 

3.2. Apertures in Finite Ground Plane 

To verify the validity of the current method when combined with GTD, a coaxial cavity with a finite 
ground plane is considered. The geometry of the cavity is shown in figure 5. The cavity is fed by a co- 
axial line with the center conductor extending up to the radiating aperture. The cavity volume is dis- 
cretized using 4474 tetrahedral elements, resulting in 4904 unknowns. For experimental verification, the 
cavity with dimensions as given in figure 5 was fabricated and its input-return loss and radiation pat- 
terns as a function of frequency measured. Figure 6 presents the input-return loss computed using the 
present method, assuming infinite ground plane, along with the experimental results. The finite ground 
plane in this case does not affect the input-return loss properties. 

The far-field patterns for the structure shown in figure 5 were calculated with equation (37) for two 
frequencies, 5.0 GHz and 9.0 GHz. Figures 7 and 8 show the calculated values along with the experi- 
mental data. The good agreement between the calculated results and the experimental data validates the 
present analysis. Because the radiation pattern of an antenna depends upon the aperture field, it is useful 
to know the tangential electric field at the aperture of the structure (shown in fig. 5). In figures 9 and 10 
the tangential electric field at the aperture is plotted for 5.0 GHz and 9.0 GHz, respectively. As 
expected, the aperture field deviates from the TEM field distribution as the frequency increases. 
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4. Conclusions 


A combined finite element method/method of moments/geometrical theory of diffraction (FEM/ 
MoM/GTD) technique was presented to analyze radiation from cavity-backed apertures with finite 
ground planes. Various numerical examples were considered to demonstrate the flexibility and validity 
of the method. The numerical data obtained using the present method and the data available in the liter- 
ature agree very well. The successful application of the GTD was demonstrated by calculating the radi- 
ation patterns of a coaxial cavity with finite ground plane and comparing the theoretical results with the 
experimental data. By virtue of the FEM, this combined FEM/MoM/GTD technique is applicable to 
arbitrarily shaped cavities filled with inhomogeneous and anisotropic materials. With proper modeling, 
microstrip antenna arrays residing in cavities with finite ground planes can be also analyzed. This tech- 
nique can also be easily extended to a finite circular ground plane. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
September 6, 1995 
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Appendix 


Modal Expansions of Feed Structures 

For completeness, expressions for the modal expansions of electric fields (eq. (19)) at the input 
plane of feed structures arc given below, with the incident field given by equation (20). 

A.l Rectangular Waveguide 

The cross section of the input plane of the rectangular waveguide is shown in figure Al. The inci- 
dent electric field is assumed to be the dominant TE 10 mode, which is given by 



(Al) 


The reflected electric fields for TE and TM modes are written with mode m = (p, q) as (ref. 10) 
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(p and q = 0, 1,2, . . . ; if p = 0, q ^ 0; if q = 0, p ^ 0) (A2) 


and 


where 


TM 

e m = xA pq { 


f pn\ [pnf flYl • \~qnf bA 

lT> os LTl x + 2jJ sm Lw + 2] 


+ y a 


(<vp 
l b , 


\ . 

\pn 

r , a\\ 


~qn 

r b\ 

sin 

[t\ 

2 + 2jJ 

cos 

T 

y + 2j\ 


+ z 


pq 




a J 


\2n 


b ) 


sin 


'pnf a 
. a V 2 


a p, = 


£ £ 
op oq 


ab 


?H2 

(p = 0, 1,2,3, ...;q = 0, 1,2,3, ...) (A3) 


1 


£7tV 


(A4) 


+ (m 

a J \b J 


and 


8 *P = 


(for p ± 0) 
(for p = 0) 


e -/ = 


(for q * 0) 
(for q - 0) 


13 



Y 


m 




> 



(A5) 


A.2 Circular Waveguide 

The cross section of the input plane of the circular waveguide is shown in figure A2. The incident 
electric field is assumed to be the dominant TE^ mode, which is given by 


e inc = A lliP 


~ J l( k p P) 00 ^ 

r 


+ §[k' J\ (k' p) sin(f)] 


(A6) 


and 



(A7) 


X 1 1 

where k' = and Y , represents the first zero of the derivative of the Bessel function 

P r c 

Equation (A9) defines Ajj with p = 1 and q-\. The reflected electric fields for TE and TM modes with 
m = {p, q ) are given by 


TE 

e m = A Pq \ P 


L P 


Z J q ( k p P)cos(q^) 


+ (})[ Ap J' q (k p p ) sin (^(f)) ] 


(A8) 


e ,™ = B Pq \^ k p J 'q (^pplsin (#)] + $ 


zJq( k p P)cos(q§) 


L P 


+ + Ym)^(^pP) sin (^)] 

* m 


(A9) 


X X 

where k' = — — and k n = — — ■ %' is the nth zero (p = 1,2, ...) of the derivative of the Bessel func- 
f r c f r c Pi 

tion J q of the first kind of order q (q - 0, 1,2, . . .). % pcj is the /;th zero {p - 1, 2, 3, . . .) of the Bessel func- 
tion J of the first kind of order q (q = 0. 1,2, . . .), and 


A - 1 

pq r ' 2 2, J ( y' ) 

^(Xpq-q) < 7 V/ W 


(A10) 


R = 1^22 1 

pq v 71 x pq J ' q (x pq ) 


(All) 
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where k = k' p for TE modes and K = kp for TM modes. 

A.3 Coaxial Line 

The cross section of the input plane of the coaxial line is shown in figure A3. The incident electric 
field is assumed to be the dominant TEM mode and is given by 


e- = p- 

inc r 


\ 


2n In 




Tine = J k o 


(A13) 


(A14) 


In a practical feed-line arrangement, the coaxial line is relatively small in diameter and the ampli- 
tudes of the higher order coaxial modes that are excited are small. A reasonable approximation is to 
assume that only a reflected TEM mode exists in the coaxial line. Hence, in equation (19), 


TE 


= P- 


2n In 


v r iy 


TM A 

e,„ = 0 


(A15) 


(A16) 


y,n = A 


(A17) 
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Normalized admittance 



▲ 

Input 


(a) Rectangular waveguide with dielectric plug radiating into half space. 



Frequency, GHz 


(b) Variation of normalized admittance with frequency for case a = 4.755 cm, b = 2.2 1 5 cm. 
L = 1.9 cm, and e r = 3.76. 

Figure 2. Rectangular aperture in infinite ground plane. 
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Norn 



Figure 3. Circular aperture in infinite ground plane. 
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Normalized admittance 



A 


Input 

(a) Coaxial line with inner radius rj and outer radius r 2 radiating into half space. 



Figure 4. Coaxial aperture excited by coaxial line. 
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Normalized radiated power, dB Normalized radiated power, dB 


FEM/MoM/GTD with finite ground plane 




Figure 7. Far- field patterns at 5.0 GFlz. 
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Normalized radiated power, dB Normalized radiated power, dB 



0, deg 

(a) Radiation pattern of coaxial cavity at cf) = 0° and cf) = 180°. 



0, deg 

(b) Radiation pattern of coaxial cavity at cf) = 90° and <f = 270°. 
Figure 8. Far-field patterns at 9.0 GFIz. 
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Figure 9. Vector plot of electric field at coaxial-cavity aperture at 5.0 GFIz. Amplitudes have been 
equalized for clarity' of presentation. 



Figure 10. Vector plot of electric field at coaxial-cavity aperture at 9.0 GFIz. Amplitudes have been 
equalized for clarity' of presentation. 
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